Light emission and finite frequency shot noise in molecular junctions: from tunneling 

to contact 
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Scanning tunneling microscope induced light emission from an atomic or molecular junction has 
been probed from the tunneling to contact regime in recent experiments. There, the intensity of the 
light emission shows strong correlation with the current /charge fluctuations at optical frequencies. 
We show that this is consistent with the established theory in the tunneling regime, by writing the 
finite- frequency shot noise as a sum of inelastic transitions between different electronic states. Based 
on this, we develop a practical scheme to perform calculations on realistic structures using Green's 
functions. The photon emission yields obtained re-produce the essential feature of the experiments. 
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I. INTRODUCTION 



When a scanning tunneling microscope (STM) tip is 
brought towards a metal surface, strong localized plas- 
mon modes develop between the tip and surface, in addi- 
tion to the propagating surface mode at the metal inter- 
face. Under an electric field, the plasmon modes interact 
with the electrons traversing the gap. This provides an 
efficient way to excite the plasmon modes electrically, and 
has become an important topic bridging nanoelectronics 
and plasmonics^KU Radiative damping of the excited 
plasmons results in light emission, which can be detected 
experimentally i n the f ar field at the same or opposite side 
of the STM tipPtiMOJ. Analyzing the emitted light can 
provide information about the nanogap; the dependence 
of light emission on the type of metal, the shape of tip 
and surface, and on the inserted molecular layer between 
tip and surface, have all been explore d 14 " 1 9 !. D ifferent 
types of plasmon modes have been detected 20 * 21 *. Most 
of these experiments are done in the tunneling regime, 
where the coupling between STM tip and metal surface 
is weak. Theoretically, it has been established that the 
excitation of plasmon modes is due to the inelastic elec- 
tronic transitions taken place near the gap 22 * 23 *. 

Recently, STM-induced light emission has been probed 
during the transition from the tunneling to the contact 
regime, both for single atom contacts and a Ceo molecular 
junctiorP^HSD. The experimental results reveal a strong 
correlation between the light emission intensity and the 
current /charge fluctuations at optical frequencies, and 
furthermore show the possibility of controlling light emis- 
sion by engineering the electronic structure. These re- 
sults at the strong coupling, contact regime seem to con- 
tradict the established theor^MH j n ^he wea k coupling, 
tunneling regime. 

A detailed modeling of these experiments needs to take 
into account the plasmon field distribution near the STM 
tip, the nonequilibrium electronic structure at high bias, 



the coupling of the plasmonic field with electri cal cur- 
rent, and the propagation of light to the far fiel d 22 * 23 ^. 
In this paper, instead of developing a full theory, we fo- 
cus on the electronic part of the problem. In particular, 
we study how the change of the electronic structure influ- 
ences the efficiency of plasmonic excitation. To this end, 
we derive a Fermi-golden-rule like expression for the finite 
frequency shot noise, and relate it to the theory of STM- 
induced light emission in the tunneling regime. We then 
express the result in terms of nonequilibrium Green's 
functions (NEGF), and develop a practical scheme to 
perform calculations on realistic structures, using infor- 
mation available from Density Functional Theory based 
NEGF (DFT-NEGF) transport calculations. We demon- 
strate how this scheme manage to capture the essential 
feature of the atomic metal and molecular contact exper- 
iments. 



II. THEORY 

In this section, we briefly summarise the the ory of 
STM-induced light emission in the tunneling regime^M] # 
Then, following Ref. [28j we introduce an approach to ex- 
press the finite frequency shot noise in a coherent con- 
ductor as a sum of inelastic electronic transitions. We 
demonstrate how the shot-noise explanation of the light 
emission in a molecular contact is consistent with the 
theory in the tunneling regime. 



A. 



Inelastic transition due to electron-plasmon 
interaction 



Following the theory of light emission from STIN/P^ESI 
and point contacts^, the interaction of the electrical cur- 
rent with the plasmon field in the tip-surface cavity is 



described by the following Hamiltonian, 
1 
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j(r)A(r)d 3 r , 
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where j(r) is the electron current density operator at 
position r. The plasmon mode, with frequency, Q, and 
spatial distribution, £(r), is represented by a vector po- 
tential, 
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Here a(a)) is the annihilation (creation) operator of the 
plasmon mode, c is the speed of the light, h the reduced 
Planck constant, and V the normalization volume. In 
principle, we may calculate the plasmon mode frequency 
and field distribution for a given a tip-surface distance. 
However, this is a daunting task for atomistic first prin- 
ciples theory and we do not consider this problem here. 
Instead, we focus only on the source of the light emis- 
sion, and investigate the effect of the non-equilibrium 
electronic structure on the emission rate. We ignore the 
spatial distribution of the mode in the x?/-plane trans- 
verse to the current, £(r) = £(z), and perform the inte- 
gration over these directions in Eq. (fTl and get 
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I(z)A(z)dz, 



= M(a + a ] ), 



(3) 



where I(z) is the surface current evaluated at z, inte- 
grated over the transverse surface. The emitted power 
from the junction is proportional to the inelastic transi- 
tion probability due to the interaction between initial (^) 
and final (^z) states originating from the tip (t) or surface 
(s) electrode, 



hf 

xn F (si - /i»)(l - n F (ef - /J,f))dsids f . 



(4) 



We employ the normalization, (i/ji\ijjj) = 5ij5(si—Sj), and 
filling given by the Fermi-Dirac distribution n F . Finally, 
we assume that the "diagonal" contributions in the z 
direction capture the main dependence of the emitted 
power on the electronic structure of the junction. Thus 
we get, 



ij 



Xnp(Si — Mi)(l — n F^f — jJLf))d£idSf. 



(5) 



This "diagonal" assumption can clearly not be justified 
per se without concrete knowledge about the spatial dis- 
tribution of the mode along with the local current oper- 
ator. On the other hand, below we will use a first prin- 
ciples method in order to calculate without any fitting 
parameters the light emission using this approximation 
and compare with the experimental trends. 



B. Current, charge fluctuations and emission rate 

Now we show that the Fermi's golden-rule rate in 
Eq. ([5| is closely related to the finite frequency shot noise 
of the electrical current, which is defined as 

«i,(o)i z ,(*)» = ((W)-(W)))(iAt)-{iAt)))), (6) 

where I(t) = e tHt / h j e -' lHt / h i s the surface current opera- 
tor along z in the Heisenberg representation and zj z' are 
two positions along the transport direction. The positive 
direction of I z is defined to be from the surface electrode 
towards the tip. Since we are dealing with the time de- 
pendence explicitly, we put the position variables z, z' as 
the sub-indices. The Fourier transform of Eq. ^ gives 
the noise spectrum, 



S Z z'(u) 



+oo 



((I z (0)Iz'(t)))e Mt dt. 



(7) 



Following Ref. [28l inserting a complete set of eigenstates 
into Eq. (m), and doing the Fourier transform, we obtain 
a golden-rule-type expression for the current noise, 



S zz ,(u>) = 27rh^JJ (^\Wf)(^f\^mS(e l -e f -fiw) 
xn F (ei ~ /ii)(l - n F (sf - /jLf))dsidsf . (8) 



The initial and final states are summed over scattering 
states from both electrodes. Equation ([8| includes both 
the Nyquist- Johnson (thermal) and shot noise contribu- 
tions. Since the energy of the emitted light is much larger 
than the thermal energy (ftw ^> fc^T), only the zero- 
temperature limit is considered. In this case, besides the 
zero-point fluctuations, the only contribution is the shot 
noise, 



S zz >(lj) = 2ttH 
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(ipt\Iz\ips}(ips\I z '\i't)de t , 



(9) 



with e s = e t —huj for positive sample bias V = V S —V t > 0. 
The "diagonal" correlation S zz gives the sum of the 
transition rates between the initial filled tip scattering 
states ipt: and the final empty surface scattering states 
ijj s , with energies e t and e s , respectively. This illustrates 
how the finite frequency shot noise can be viewed as 
inelastic electronic transitions between the tip and sur- 
face scattering states. Equation. (J9| is closely related 
to Eqs. ppl), and consistent with the theory for STM- 
induced light emission in the tunneling regime. The pos- 
itive frequency/energy part of the noise spectrum corre- 
sponds to the photon emission, relevant to the experi- 
ment, and the negative part to the absorption process. 
We notice that if z and z' are located at the surface and 
tip electrode, respectively, then according to charge con- 
servation, 



Id = Qd = I z - I z ' 



(10) 



the charge fluctuation in the central molecule/" device" 
region(d), 



&dd — &zz H" &z' z' &zz' &z 



(ii) 



Similarly for the fluctuation of the average current I a = 

^aa = T \^zz H" ^z'z' + ^>zz' H" ^z'z) • (12) 



III. NUMERICAL SCHEME 

We aim at a formulation targeting the DFT-NEGF 
approach to atomistic electron transport, such as the 
SIESTA/TranSIESTA method 30 and similar methods 
employing a localized basis set. In these the whole system 
is separated into a central device region(d), and two elec- 
trode regions, here the tip (t) and surface (s) electrodes. 
The electrodes are represented by the self-energies. In or- 
der to directly employ the DFT-NEGF formalism we will 
rewrite Eq. ([9| in terms of the device Green's functions 
and the self-energies (E s ,E t ) folded into the same device 
region representing the coupling of the device region to 
tip and surface electrodes, respectively. By our choice 
of device region we effectively define separating surfaces 
between the regions. 

As an example we now consider the current evaluated 
at the surface electrode. In order to calculate the surface 
electrode current fluctuations, S ss (u), an explicit expres- 
sion for the surface current is needed in terms of quanti- 
ties readily available in the DFT-NEGF calculation. The 
current matrix I S: can be written aP^, 



(V ds - V sd ), 



(13) 



where Vds is the coupling matrix between the device and 
surface electrode, Vj s is the complex conjugate of the ma- 
trix Vds, an d e is the electron charge. We assume an or- 
thogonal basis set; however, a generalisation to the non- 
orthogonal case is straightforward by a Lowdin transfor- 
mation. Next, we use the Lippmann-Schwinger equation 
to make connection between the scattering states and the 
retarded Green's functions, G(e), 



\Me)) = \Me)) + G(e)V T \Me))- 



(14) 



Here \i/j 3 (e)) and \(f> 3 (e)) are the scattering states from 
the semi-infinite surface electrode with and without cou- 
pling to the device, respectively. Note that (j) s is non-zero 
only in the surface electrode, but ip s spans over the whole 
region including both electrodes and the device. The 
coupling matrix, Vr, represent the coupling between the 
device and the two electrodes, while G(e) is the retarded 
Green's function of the whole system including the effect 
of Vr- Using the projection matrix P t + Pd + P s = I, 
and the fact that Vr and \<j) s ) are non-zero only at spe- 
cific region of the whole system, we can write the current 



matrix element as, 
(^(e)|7 8 |^(e_)> 



te 
H 



{i>t(e)\W a (s.,e)\i> d s (e-)) 



(15) 



(16) 



defining e_ = e — fku, \tps( e )) = Pdl^sis)), and 

Wi(e-,e) = G- d \e.) + £;( £ _) - sj( £ ). 

Note that all quantities are projected to the device region 
and thus depend on the actual splitting into regions. 

Finally, we can write the transition rate at zero tem- 
perature as, 

S ss {lo) = [Tr [W a {e-,e)A a (e-)W}(e-,e)A t (e)] de , 
Je 

(17) 



where the integral is defined as 

^2 



r\eV\/2 
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/ • de = 0(eV - hcu)—- / • de, (18) 

Je ZttH Jhoj-\eV\/2 

with 6{x) being the Heaviside step function, A s (e) = 
G d (e)T s (e)G d (e) = 2nJ2 i=s \^{e)){^{e)\ is the device 
spectral function due to scattering states from the surface 
electrode, similarly for A t , and T s = z(E s — E|). In the 
same way, we get the fluctuations of the tip current, 



St, 



(w)= / Tr\w}(e,e-)A s (e-)Wt(e,e-)A t (e) 
Jo L 



ds , 

(19) 
and their cross correlation, 
S st (oj) = S; s (u) (20) 

' Tv[W s (e-,e)A s (e-)W t (e,e-)A t (e)]de . 



I, 



Equations (17 [21) are our main formal results, where we 
have written the finite frequency shot noise in terms of 
the Green's functions and self-energies, readily available 
from DFT-NEGF calculations. The difference between 



Eqs. (17) and (19) reveals the position dependence of 



finite frequency noise. Importantly, they both yield the 
standard result in the zero- frequency limit 32 . 

Assuming constant self-energies (E s ,Et), and decou- 
pled transmission eigenchannels 31 at different energies, 
we arrive at more physically transparent expressions, 

S ss (u) = V / T n (e)(l-T n (e-))de, (21) 

Sttiw) = V / T n (e-)(1-T n (e))de, (22) 
n J e 

for positive sample voltages. (See Appendix [A] for the 
full result of S ss (uj) at finite temperature). The two ex- 
pressions are exchanged for negative bias. Here T n are 
the channel transmissions calculated for a particular bias. 
Unfortunately, we are not able to write the cross corre- 
lations S st and Sts in terms of the eigen-transmissions 



IV. RESULTS 

Now we apply the method outlined above to calculate 
the light emission from the STM resembling two recent 
experiments where the tip is brought into contact with 
(i)a Ag adatom on a Ag(lll) surface^, and (ii) a Ceo 
molecule a Cu(lll) surface 26 . In the experiments, two 
type of photons with energy smaller and larger than the 
applied bias are detected. They are attributed to one- 
and two-electron process, respectively. Here, we focus 
only on the former. We used the SIESTA/TranSIESTA 
cod d 30 * 33 -! with the generalized gradient approximation 
(GGA-PBE) for exchange and correlation 3 -^. For the Ag- 
system, we use a single-^ polarized basis-set for the Ag 
atoms. For the Ceo -system, we use a double-^ basis-set 
for the carbon atoms, and a single-^ basis-set for the bulk 
electrode Cu atoms. For both systems, to accurately de- 
scribe the surface and/or the chemical bonding with the 
Ceo , an optimized diffuse basis set was applied for sur- 
face layer atoms and the tip 35 
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A. Ag adatom on Ag(lll) 

Figure [IJa) shows a subset of the structures used 
in the calculations, going from tunneling to contact 
regime. A 4 x 4 surface unit-cell were used, and we use 
2x2/5x5 surface k-points to sample electronic struc- 
ture/transmission. We relaxed the two surface layers, the 
tip and the adatom at zero bias. After the relaxation, 
transport calculations were done for a bias of V = ±1.5 
V. Figure [lib) shows the transmission eigenchannels for 
the structures in Fig. []Ja). From Fig. [TJb) it is evident 
that, (i) there is only one dominate transmission eigen- 
channel, and (ii) there is a small asymmetry in the trans- 
mission for the two bias polarities. Figure lllc) shows the 
change of the average conductance when going from tun- 
neling to contact on a log-scale. In the tunneling regime, 
the conductance depends exponentially on the tip-atom 
distance, while it develops to a plateau upon contact as 
typically seen in experiments 25 . 

The emission rates was evaluated for a plasmon energy 
of h£t = 1.2 eV using Eq. ([5|, or equivalently Eqs. (17- 
19). In order to map out the spatial distribution, the 
emission rate is calculated for the surface current de- 
fined at 6 different interfaces, shown in Fig. [2Fa)-(b). 
From these calculations, we observe that the emission 
rate does not change significantly for interfaces in the 
same electrode, while they are quite different for the two 
electrodes, and for the tip- adatom interface. 

To relate the emission rates to the intensity of light 
emission, we need to do an average of the surface cur- 
rents, taking into account the spatial distribution of the 
plasmon mode, £(z). Since we do not have specific knowl- 
edge about the mode we will choose to do it in the sim- 
plest possible way here. Firstly, we take the equally- 
weighted average of all the surface layers. In reality, 
the mode distribution may change with the tip-surface 



FIG. 1: (a) A subset of structures used in the calculation, 
going from tunneling to contact. In the final structure, one 
tip atom is pushed aside when forming contact. The two 
surface layers, the tip and the adatom are relaxed at zero bias 
for each structure. The numbers show the distance between 
the two fixed layers and between the tip-adatom in unit of 
A. (b) Transmission eigenchannels at V = V s — Vt — ±1.5 
V, going from tunneling to contact (top to bottom), for the 
structures shown in (a), (c) The average conductance as a 
function of surface layer separation, showing the transition 
from tunneling to contact. 



distance. In the tunneling regime, we expect a high 
weighting-factor in the region between the tip-surface 
gap. On the other hand, upon contact, due to the high 
conductance, we expect the fi eld distribution to spread 
out into both electrode d 36 * 37 *. Secondly, as mentioned 
above we will use Eq. (I5| instead of Eq. Q, so we ig- 
nore the cross terms involving surface current at differ- 
ent positions. The final results for the two bias polari- 
ties are shown in Fig. [2jc)-(d) together with the approxi- 
mate calculation using Eqs. ( |21p2 ), and the result using 
zero- frequency noise which was employed in Ref. [25l We 
see that the qualitatively trends are similar for all these 
calculations. We see a plateau in the tunneling regime, 
and the development of a dip at contact around the fully 
transmitting single channel for G = IGo, consistent with 
the experiments 25 . 



B. C 60 on Cu(lll) 

For the Ceo system, a 4 x 4 surface unit-cell, and 
2 x 2/10 x 10 surface k-points was used in order to sam- 
ple the electronic structure/transmission. Here due to 
the surface reconstruction in the experiment J 26 * 38 * the two 
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FIG. 3: (a) All structures considered in the calculation. In the 
6th a deformation of the tip occurred. It is disregarded in the 
following calculations. The two surface layers, Ceo and the tip 
were relaxed at zero bias for each electrode separation, (b) 
Transmission eigenchannels at V = ±1.5 V for the structures 
shown above. 



FIG. 2: (a)-(b): Calculated emission rates from Eq. ^ for 
I z defined through 6 different surfaces, shown above, for plas- 
mon energy Ml = 1.2 eV. (c)-(d): Calculated yields from 
average of the emission rates at different surfaces (squares). 
Also shown are the results from av erage o f S ss and S t t using 
the approximated expressions Eqs. ( 21|22 ) (circles), and from 
the zero- frequency noise calculation used in Ref. 25 'triangles). 
All of them give qualitatively similar results. 



first surface layers and tip were relaxed at zero bias to 
0.02 eV/A at different tip positions. Thus, we do not cap- 
ture the abrupt jump-to-contact observed in the experi- 
ment at finite negative bias in our calculations. Figure [3] 
shows the five different structures considered in the cal- 
culations, together with the transmission eigenchannels 
at V = ±1.5 V. Different from the Ag system, when mak- 
ing the contact, there are now mainly three contributing 
eigenchannels. 



Importantly, the emission rates for the Ceo systems are 
quite different for the two bias polarities (Fig. EFa)-(b)). 
For positive sample bias, the magnitude at 4 different 
surfaces is comparable. But for the negative bias, the 
fluctuations near the surface electrode are 4 times larger 
than that of the tip electrode. This can be explained 
as a consequence of the appearance of the HOMO level 
in the bias window, as discussed in Ref. [26j When this 
happens, the charge in the HOMO state begins to fluctu- 
ate. This generates new available final states for inelastic 
transitions, which contribute to high frequency noise at 
the plasmon frequency. Since the molecule couples better 
to the surface than the tip, the charge fluctuations are 
compensated mainly by the surface- current fluctuations. 
Consequently, the calculated yields show different trends 
at negative and positive bias when going from tunneling 
to contact, as shown in Fig.[4^c)-(d). Notice that the ap- 
proximate results do not capture this polarity dependent 
yields. Neither does the zero- frequency noise calculation. 
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FIG. 4: (a)-(b) Similar to Fig. [2ja)-(b), calculated emission 
rates at 4 different surfaces for the Ceo system using hQ = 1.2 
eV at V = ±1.5 V. (c)-(d) Similar to Fig. |2] (c)-(d). 



V. CONCLUSIONS 

We have developed a practical scheme to calculate 
the finite-frequency shot noise of the electrical current 
through a coherent molecular conductor within a DFT- 
NEGF approach. By a spatial average, we re-produce 
qualitatively the essential features of two recent experi- 
ments, confirming the hypothesis that the current/charge 
fluctuations are the energy source of STM-induced light 
emission from molecular junctions, going from tunnel- 
ing to contact. Furthermore, by writing the shot noise 
expression into a Fermi-golden-rule form, we have estab- 
lished a connection with the theory of light emission in 
the tunneling regime, based on inelastic electronic tran- 
sitions. 



Here, we have focused on the source of the light emis- 
sion, which is the inelastic electronic transitions induced 
by current. However, to get a quantitative understand- 
ing of the experimental results, in a semi-classical model 
of the electron-plasmon coupling, the following questions 
have to be addressed: (1) the spatial field distribution 
of different plasmon modes near the STM tip, (2) their 
detailed coupling with the current. These questions are 
also important if we want to distinguish the localized gap 
mode from the propagating surface mode. Recent exper- 
iments showed that the tunneling electrons can couple 
to both types. An alternative way to proceed is to per- 
form time dependent DFT calculations. So far, model 
structures have been considered^! with this approach. 
However, it is very challenging to perform calculations 
on realistic structures involving coupling to the metallic 
surfaces in order to approach the experiments. 



Appendix A: Frequency dependent noise at finite 
temperature 

At finite temperature, to evaluate the surface current 
correlation, we need all the matrix elements. The other 
three read 



(Me)\I s \Me-)) 
(Mz)\Is\Me-)) 



% -%(Me)\W!(e,e-)\Me-)), 
-(V> t (e)|E.(e-)-Ej(e)hMe-)>, 



(Me)\Is\Mz-)) = ¥ (^( £ )|Sl(e)-S t (e_)- w /|^( £ _)). 

Assuming a constant self-energy, for positive sample bias, 
we have the full result for surface current noise at finite 
temperature 

e 2 

a(3 

with 

Ctt{u>) = JTr[T{e)T{e-)]Anfde, 

C ss (uj) = I Tt[(ujI -ir t )A a (e-)(uI + iT t )A s (e)] An s F s de, 

Cst = J Tr[(I -T(e))T(e-)]An^ds, 

Cts = J Tr[(I - T(e_))T(e)] Ar#cfe, 



where 



a(3 
l F 



An£ p = n F (e,iL a )(l -n F {e-,np)). 



The above result includes both the Nyquist- Johnson 
(thermal) and the shot noise. Notice the different form of 
C ss from Cu- It is related to the complex reflection coeffi- 
cients in the scattering approach discussed by Biittikei 40 . 
Physically, it means that even when the transmission is 
zero, there still could be fluctuations at the surface elec- 
trode at finite temperature. 
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